Possibility of s-wave pion condensates in neutron stars revisited 
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We examine possibilities of pion condensation with zero momentum (s-wave condensation) in neutron stars 
by using the pion-nucleus optical potential U and the relativistic mean field (RMF) models. We use low-density 
phenomenological optical potentials parameterized to fit deeply bound pionic atoms or pion-nucleus elastic 
scatterings. Proton fraction (Yp) and electron chemical potential (/le) in neutron star matter are evaluated in 
RMF models. We find that the s-wave pion condensation hardly takes place in neutron stars and especially has 
no chance if hyperons appear in neutron star matter and/or bi parameter in U has density dependence. 

PACS numbers: 21.65.Jk, 26.60.-c, 36.10.Gv, 25.80.Dj, 



Pion condensation in neutron stars has a long history of 
study from the first suggestion in early 1970s by Migdal iQl] 
and Sawyer The pion-nucleon interaction is attractive in 
p-wave, then the main interest along this line was to explore 
possible appearance of pionic excitations with zero energy and 
finite momentum, i.e. p-wave pion condensation in nuclear 
matter. Possibilities of p-wave pion condensation in finite nu- 
clei had been investigated extensively in 1970's and 1980's. 
Those possibilities were denied by the non-observation of 
anomalous angular momentum distribution in the inelastic ex- 
citation of the pionic quantum numbers [4]. Possibilities of p- 
wave pion condensation at high densities were also considered 
to be improbable based on the universal repulsion assumption, 
9'na 9'nn ^ ^-^ ' ^-^ (^fl- 1990's, new experiments on 
the Gamow-Teller giant resonances were performed, and the 
sum rule value including the 2p — 2h states was found to be 
around 90 % i^, suggesting that the transition to the A re- 
gion is weak and pJvA ^ '^■^ would be smaller than gJvA? 
In addition, microscopic variational calculation fsf] suggests 
7r° condensation in symmetric nuclear matter at high densi- 
ties (pb > 0.2 fm^"^) generated from the A mediated three 
nucleon force. Thus at present, we cannot completely deny 
the possibility of pion condensation in dense matter, and it is 
necessary to examine all the ingredients of ttN and 7r-nucleus 
interactions with updated experimental and theoretical knowl- 
edge. 

The study of the in-medium pion properties, especially the 
s-wave pion-nucleus interaction, has been recently developed 
in experiment. Precise observations of deeply bound atomic 
states of TT" in Pb and Sn isotopes jgl ll^JjU [l2i [l3ll and low 
energy pion-nucleus elastic scattering ||14ll provide us with 
detailed information of density dependent optical potentials 
at low densities. Theoretical calculations of in-medium pion 
self-energy have also experienced much progress based on 
chiral dynamics |[ll[ll[ll[ll- 

Together with these developments, it may be interesting to 
revisit pion condensation at high densities, such as in neutron 
stars. Motivated by the recent progress in the s-wave pion- 
nucleus interaction, we concentrate on the study of pion con- 
densation with zero momentum (s-wave condensation), for 
simplicity. Even such a limited investigation would make 
progress of our understanding of dense matter. The s-wave 



pion condensation takes place in neutron stars with nuclear 
matter instability where the transition n ^ pn^ becomes en- 
ergetically possible lfl9ll . This happens with — /ip > E-^, 
where and jip are the neutron and proton chemical poten- 
tials, respectively, and is the 7r~ energy at rest. Due to (3 
equilibrium under the charge neutral and neutrino-less condi- 
tions, /i„ — /ip can be written as the electron chemical potential 
/ie = Pn — Mp- Since nn interaction is more repulsive than 
pn, p„ is pushed up then /ie increases compared to Fermi gas 
value in neutron rich matter. For example, relativistic mean 
field (RMF) models suggest |[l9[ I20I1 that /Xe largely exceeds 
the in-vacuum tt^ mass at nuclear density Pb = {1 ^ 5)po in 
neutron stars, where po is the saturation density. 

In this paper, we examine whether the condition for the 
s-wave TT^ condensation, E^^ = Pe, is satisfied in neu- 
tron stars for E^ and p^ obtained in our present knowl- 
edge of low density pion optical potentials and equation of 
state (EOS) of neutron star matter. For Et^, we use various 
pion optic al p otentials fitted so as to reproduce the pionic 
atom lll2i I211 122! l23i l24l I25I1 and 7r-nucleus elastic scatter- 
ing data II14I1 . Although the fitted optical potentials in nor- 
mal nuclei may not be extrapolated to higher density and/or 
highly asymmetric nuclear matter, it is interesting to exam- 
ine the present status and to think about next steps. For pe, 
we adopt the results calculated in RMF with several param- 
eter sets ll26l I27I HH l29tl . which explain the bulk properties 
of nuclei such as the binding energy and the charge radius in 
a wide mass range. We also use proton fractions (1^) evalu- 
ated with RMF for calculating pion optical potentials' . There 
are several theoretical works on the s-wave pion condensation 
at higher densities ll30ll . while the connection to low density 
phenomena observed in experiments is not clear yet. Since 
we concentrate on the possibility of the s-wave tt" conden- 
sate, we do not consider the double pole condition for 7r~7r+ 
pair creation, which takes place with finite momentum. 

The TT^ energy in uniform matter may be evaluated by 
Etj = \Jm\ + 2mT^U (p = 0) with the real part of the 



A preliminary study has been done along the same line for limited combi- 
nations of pion potentials and an RMF parameter set I20I1 . 
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TABLE I: Pion potential parameters. The upper four sets by the pi- 
onic atom are taken from Ref. IToll . in which 60 and BJ™ were read- 
justed to reproduce the recent data of the deeply bound tt" states in 
Pb with fixing the other parameters as the ori gina l values given in 
Ref. dH] for T (H for BFG, IH for SM and El for ET. 





parameter 


bo 


61 


T-iRe 

-Bo 


a 


system 


set 


(JUtt ) 


(rriTT ) 


yrriTT ) 






T 


-0.034 


-0.078 








pionic 


BFG 


-0.025 


-0.085 


-0.021 





atom 


SM 


-0.027 


-0.12 










ET 


-0.020 


-0.0873 


-0.049 







NOG [25J 


-0.013 


-0.105 










KYQ2] 


-0.0233 


-0.1473 


-0.019 


0.367 


pion-nucleus F-Cri41 


-0.009 


-0.114 


-0.040 





scattering 


F-WQ4] 


-0.009 


-0.081 


-0.040 


0.391 




Asym. Nucl. Matter (¥^=0.2) 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 
Pb (fin"') Pb (fm"') 



E, T 
BFG 
SM 



ET 

NOG 

KY 



F-C 

F-W - - 
MOW 



KKW - 



FIG. 1: (Color online) Et^ in symmetric (Yp = 0.5, upper) and asym- 
metric (Yp = 0.2, lower) nuclear matter. 



ing self-energy, abbreviated as MOW 13, 



energy-independent potential U based on Ericson-Ericson pa- 
rameterization Olll: 



U = -- 



2tt 



(1 + e) (6oP.+5i5p)+(l + |)s«>(2) 



with e = ttl^/Mn, Pb ^ Pn + Pp, 5p = /?„ - pp = Pb{1 - 
2Yp) and a squared density p*^^) defined below. This potential 
is related to the pion self-energy via = 2to^(J7 + iW) 
with an imaginary potential W. The s-wave ttN potential pa- 
rameters (bo, bi,B^) was determined with special care from 
precise measurements of pionic atom and pion-nucleus scat- 
tering data. 

In Table m we summarize the parameter sets adopted here. 
The upper 6 sets were determined from the pionic atom data 
and the lower two from pion-nucleus scattering. For the for- 
mer parameter sets except NOG, &o = ^0 p'^' = p| are 
used, whereas, for the latter, double scattering modifications 
were explicitly included by bo — Bq — 3(6q + 2bi)kp / {2tt) 



with kp = (37rVB/2)^/^ and p' 



(2) 



Pe 



Sp is used. 



For NOG, 60 = 60 + Sbo - 3(1 + e){bl + 26?)fci./(27r) 
with 5bo = —0.0053™^^. In the parameter sets of KY 
and F-W, 61 is assumed to have density dependence through 
bi = bi/{l — apg/po) 113. mtl with finite a and po = 
0.17 fm-^ This is a consequence of the pion wave function 
renormalization associated with energy dependence of the op- 
tical potential UtIIssIi and the renormalization was performed 
at = TOtt. For more realistic calculations in dense nuclear 
matter, the wave function renormalization should be done at 
the in-medium pion mass and other parameters should be also 
renormalized. 

The phenomenological potentials are determined with a 
fixed pion energy. The energy dependence of the optical 
potential can be estimated by theoretical calculations. The 
s-wave in-medium pion self-energy ^.^^{Et^), equivalent to 
the optical potential, was derived based on the chiral per- 
turbation theory in Ref. |fl5|] within a linear density approx- 
imation. The in-medium pion mass is obtained by solving 
TO*-^ — TO J — E(to* ) = 0, which automatically takes account 
of the wave function renormalization. Here we use the follow- 



2pE 
P 



C2 + C3 



' 2/2 



with ci = -0.81, C2 = 3.20, cg = -4.66 in units of GeV^i 
and / = 88 MeV. We also consider the tt^ self-energy cal- 
culated by an in-medium chiral perturbation theory in 0{p^) 
discussed in Ref. [[TtIi (KKW), which reproduces well the en- 
ergies and widths of deeply bound tt^ atomic states in Pb. 

Let us see the model dependence of the pion energy E.^ . In 
Fig.[Tl we show i?^ in the cases of symmetric (Yp = 0.5) and 
asymmetric {Yp = 0.2) nuclear matter. The proton fraction 
= 0.2 is a typical value in neutron star matter obtained in 
RMF, as shown later in Fig.|2l The negative signs of the coef- 
ficients (bo, bi , B^°) imply that tt^ feels repulsive potential in 
nuclear matter The phenomenological pion potentials in sym- 
metric nuclear matter agree well with each other at low den- 
sities below Po. whereas, in asymmetric nuclear matter with 
Yp ~ 0.2, we have 50-100 MeV ambiguities at po- In 

order to fix the large ambiguity of the potentials in asymmet- 
ric nuclear matter, it is very interesting to obtain pionic atom 
and scattering data in neutron rich nuclei 134]] . 

RMF models have been developed to describe bulk prop- 
erties of nuclei and nuclear matter with the mean field via 
the meson fields. We here adop t the RMF models, NLl ll26ll . 
NL3 TMl iH, SCL having the Lagrangian in the 
following form, 



{\g3(T^ + \gia^ (NLl, NL3, TMl) 
1 -a./scL(^//.) (SCL) 



where ip, a, oj, p represent nucleon and a-, ui- and p-meson 
fields, respectively, and /scL(a^) = log(l— x)+x+x2/2. The 
model parameters are summarized in Table |ll] These RMF 
models describe the binding energies of heavy semi-double 
magic nuclei well, and are expected to give reasonable EOS 
of nuclear matter. We have solved the (3 equilibrium condition 
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TABLE II: RMF parameters. In SCL, gs, and are from the expansion of /scl- 





g<TN 


guN 


9pN 


53(MeV) 


54 




m<,(MeV) 


m^(MeV) 


mp(MeV) 


NL1[2&] 


10.138 


13.285 


4.976 


2401.9 


-36.265 





492.25 


795.359 


763 


NL3[2T| 


10.217 


12.868 


4.474 


2058.35 


-28.885 





508.194 


782.501 


763 


TM1[28J 


10.0289 


12.6139 4.6322 


1426.466 


0.6183 


71.3075 


511.198 


783 


770 


SCL|2a 


10.08 


13.02 


4.40 


1255.88 


13.504 


200 


502.63 


783 


770 



0.3 

a 0.2 

0.1 

^ 
> 

I 200 

i 

W 100 



Neutron star matter. . . 


NLl 

NL3 

TMl 

SCL 

TMl-SM 

lOTSY 




E/B 



0.2 0.4 J 0.6 0.2 0.4 ^ 0.6 0.8 

Pb (fm" ) Pb (f™ ) 



FIG. 2: (Color online) RMF results of proton fraction (V^ = pp/ pB, 
upper panel), electron chemical potential {p,e, lower left panel) and 
energy per baryon (E/B, lower right panel). 



in cold neutron star matter. 



Me 



Pp 



M 



N 



kp + 9oj^ T gpP : 



(1) 
(2) 



where M 



N 



M 



N 



gacr represents the effective mass of 
nucleon. As shown in Fig. |2] calculated values of E/B, Yp 
and /ie in neutron star matter are consistent at low densities 
(Pb < Pq), since meson-baryon coupling constants are well 
determined by the binding energies of heavy-nuclei. Sig- 
nificant differences are found in E/B at higher densities, 
where the mesons have large expectation values and the self- 
interaction terms {(t^ , a'^ , oj'^) contribute to E/B consider- 
ably. While we have small differences in Yp and pe, the model 
dependence is smaller compared with those in E/B and E^. 
As we can see from Eq.(2), pn — Pp is modified from the 
Fermi gas value with M* by the p meson, whose coupling 
with nucleons is well constrained by nuclear binding energies, 
and higher order terms of the p meson are not included in the 
RMF models under consideration. As a result, model depen- 
dence of the isospin dependent potential, gpp, is only around 
10 MeV at ps = 0.8 fm~^. In Fig.|2]we also show the results 
of some RMF models including hyperons (TMl-SM Issll and 
lOTSY lEoll ). With hyperons, the proton fraction and electron 
chemical potential significantly decrease. 

Now let us compare the electron chemical potential p^ 
and the pion energy i?^ as functions of pg (Fig. O. The 
RMF results of TMl are adopted for the proton fraction 
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FIG. 3: (Color online) Electron chemical potentials in RMF models 
and pion energy in neutron stars. 



(1^) to evaluate Ej^. Results with hyperons denoted by 
lOTSY are also shown. The left panel of Fig. [3] shows the 
comparison of pe and i?^ obtained from the pionic atom 
data, namely, Tauscher (T), Batty-Friedman-Gal (BFG), Seki- 
Masutani (SM), Ericson-Tauscher (ET) and Kienle-Yamazaki 
(KY). In these potentials, density dependence of the poten- 
tial parameters is not taken into account except for KY. The 
right panel of Fig. [3] compares p^ with E-,^ obtained from the 
pion-nucleus scattering data and the theoretical calculations. 

We find that £^7r obtained from the potentials with density- 
independent bi are very close to pe at high densities pg > 
0.3 fm^'^, and in further dense nuclear matter, p^ exceeds 
i?7r obtained with some of the parameter sets. We could have 
possibility for the s-wave pion condensation to take place in 
dense neutron star matter. However, since inclusion of hyper- 
ons makes pe suppressed, E-^ is found to be larger than pe 
(lOTSY) in most cases. Thus the s-wave pion condensation 
would not take place, if hyperons could participate in neutron 
star matter. 

We note that even more attractive hyperon potentials make 
Pe smaller (TMl-SM). Also in non-relativistic variational 
treatments |H, [36ll . symmetry energy and pe are generally 
smaller at Pb > po than in relativistic models. In the case 
of the density-dependent bi, which is a consequence of the 
renormalization of the pion wave function and is required to 
explain the pionic atom data of Sn isotopes |fT2f|, the pion 
self-energies are more repulsive. Nevertheless, it is impor- 
tant to note that, as already mentioned, the renormalization of 
the wave function has been done only for the bi parameter in 
linear 6p. Thus, for more quantitative discussion, it is nec- 
essary to improve the phenomenological and theoretical pion 
optical potentials in a consistent way, for instance as done in 
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Ref. ifisll . It is also desired to include the effects of short- 
range and tensor correlations on /^jg under /3-equilibrium [8] 
in relativistic frameworks js^ . 

In summary, we have discussed the in-medium pion energy 
in the context of possibility for the s-wave pion condensa- 
tion to take place in neutron stars. We have compared the in- 
medium pion energies determined from pionic atom or pion- 
nucleus scattering data with the electron chemical potential 
evaluated in relativistic mean field (RMF) models, using the 
RMF result of the proton fraction. With our present limited 
knowledge of the in-medium pion properties obtained in ex- 
periments, we could conclude that the s-wave pion condensa- 
tion would not take place in neutron stars with hyperons. It is 
certainly necessary to investigate in-medium pion self-energy 
theoretically in more elaborated prescription to go beyond nu- 
clear density. Especially energy dependence of the pion self- 
energy should be treated in more proper ways for higher den- 



sities. At the same time, experimental observations of pionic 
atoms and scattering in neutron rich nuclei are essential to fix 
ambiguities in the pion optical potentials in asymmetric nu- 
clear matter. Precise knowledge of pion self-energies at high 
density is also important to study finite temperature process 
such as black hole formations, where the temperature can be 
as large as T = 70 MeV ifssll . At such high temperatures, pion 
contribution could be significant depending on the in-medium 
pion mass. 
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